Finite Temperature Ordering in the Three-Dimensional Gauge Glass 
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We present results of Monte Carlo simulations of the gauge glass model in three dimensions using 
exchange Monte Carlo. We show for the first time clear evidence of the vortex glass ordered phase at 
finite temperature. Using finite size scaling we obtain estimates for the correlation length exponent, 
v = 1.39 ± 0.20, the correlation function exponent, r\ — —0.47 ± 0.07, and the dynamic exponent 
z = 4.2 ± 0.6. Using our values for z and v we calculate the resistivity exponent to be s = 4.5 ± 1.1. 
Finally, we provide a plausible lower bound on the the zero-temperature stiffness exponent, 9 > 0.18. 
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I. INTRODUCTION 

High-temperature superconductors have a phase dia- 
gram that is rich with physically diverse phenomenal!]. 
In the mixed state of a pure type-II system one finds 
the Abrikosov latticea of triangularly arranged vortices. 
This vortex lattice prohibits superconductivity because 
any perpendicularly applied current produces a Lorentz 
force causing the vortices to move, dissipating energy. 
The addition of disorder, however, drastically changes 
the behavior of the mixed state, p. Correlated disorder, 
such as from heavy ion irradiationEI or twin boundaries,!!! 
causes the vortex lines to locally align with and adhere to 
the defects; this destroys the long-range order of the lat- 
tice and produces,-©, superconducting glassy phase known 
as the Bose glassE. Random point disorder, e.g., from 
proton irradiation!], pins the vortices in random positions 
creating a diffcrcntjJAfpe of superconducting phase known 
as the vortex glassQtl. 

Similar to a spin glass, the ordered state of the vortex 
glass is characterized by the phase of the superconducting 
order parameter randomly oriented in space but frozen in 
time. A system must have a Hamiltonian that contains 
both randomness and frustration in order to exhibit this 
type of behavior. One such model that has been used 
extensively |-to simulate the vortex glass transition is the 
gauge glassoo. 

The gauge glass is important to stud;y_.because like 
other spin glasses in 3 spatial dimensions Ji3~0 the exis- 
tence of a finite temperature transition within the gauge 
glass is controversial. The first estimate of the critical 
temperature was given in 1990 by Huse and SeungtH, 
< T c < 0.6. A year later thiS|-estimate was im- 
proved significantly by Reger et al.,ll3 who found that 
T c — 0.45 ±0.05, but their data was insufficient to estab- 
lish that spin glass Afder occurs below this temperature. 
Subsequent studiesotj have focused on defect energy 
scaling to determine the existence of a finite T c by cal- 
culating the stiffness exponent 9. Unfortunately, these 
results seem to be dependent on the chosen boundary 
conditions, and there is little agreement on the values of 
9. Consequently it is not clear whether the lower critical 



dimension is 3 or less than 3. 

In this paper we attempt to settle the controversy 
by studying the gauge glass using exchange -Monte 
Carlo,EJE!J also known as parallel tempering.EjCfl This 
technique allows us to simulate lower temperatures and 
larger sizes than previously possible. The main features 
of our work are as follows. 

1. We present for the first time clear evidence of 
vortex-glass ordering at finite temperature in the 
three-dimensional gauge glass model. We. present 
superior data to that previously available^!! and ob- 
tain a more accurate determination of the critical 
temperature T c /J = 0.47 ± 0.03. 

2. Using finite size scaling we determine the correla- 
tion length exponent v = 1.39 ± 0.2jQ-.ta a higher 
degree of precision then earlier works .P°H 13 I The cor- 
relation function exponent is also obtained, r\ = 
—0.47 ± 0.07. To our knowledge, this is the first 
numerical estimate of r\ for the gauge glass from 
simulations. Assuming hyperscaling, these two ex- 
ponents completely determine the universality class 
of the gauge glass model. 

3. Using standard Monte Carlo, we determine the dy- 
namical exponent z — 4.2 ± 0.6 and compare our 
results to experimental measurements of the resis- 
tivity exponent s = v{z — 1). We find s = 4.5 ± 1.1. 

The layout of the paper is as follows. In §|ll] we de- 
scribe the model while in §|l| we discuss the observables 
that we measure. In §|w] we discuss our implementation 
of the exchange Monte Carlo method and our tests for 
equilibration. Our results for statics are discussed in 
while our results for dynamics are given in § VI. In 



we summarize our results and give some perspectives for 
future work. 



II. THE MODEL 

The gauge glass describes the physics of a disordered 
type-II superconductor at distances larger than some 
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characteristic length scale A, .beyond which order in the 
flux lattice has been broken.ElO One can then imagine 
the system as a granular superconductor (in an applied 
magnetic field) with an inter-grain separation of order 
A. Such a system can be modeled as an-array of nearest 
neighbor coupled Josephson junctions.BEj This leads to 
the Hamiltonian 



TL = — J cos( 



4>j - Aij), 



(1) 



where each site i on an N = L x L x L cubic lattice has an 
associated phase angle (pi, J is a positive ferromagnetic 
(Josephson) coupling between nearest neighbors, and Ay 
is proportional to the line integral of the vector poten- 
tial along a straight line path from site i to its nearest 
neighbor site j, 



4>o 



A ■ dl. 



(2) 



$o = hc/(2e) is the flux quantum. 

The Hamiltonian of the gauge glass is given by Eq. (|l|) 
in which the Ay are quenched random variables uni- 
formly distributed from to 2ir. Note that, by contrast, 
restricting the Aij only to values and ir leads to the 
XY spin glass, which is equivalent to setting Aj = 
and taking the interactions to be ± J at random. 

The gauge glass is perhaps the simplest model of a dis- 
ordered type-II superconductor that contains the correct 
order parameter symmetry as well as the randomness and 
frustration necessary for glassy behavior. However, there 
are some features it ignores. U 

The model ignores screening since the Aij in Eq. (jl|) 
are quenched; there are no thermal fluctuations in the 
magnetic field. This corresponds to tie extreme type- 
II limit in which the Ginzburg-LandauEH parameter n — 
A/£ — > 0Oj where A is the penetration depth and £ is the 
coherence length. This limit may be realistic since k can 
be quite large in high T c superconductors, e.g., k rs 90 
for YBa2Cu307_a. It seems, however, that when the 
gauge glass is modified to include strong screening, the 
finite temperature transition to the vortex glass phase is 
rounded putji three dimensions very close to the pu- 
tative T c r3~E3 This rounding probably takes place over 
such a narrow temperature region that will be very diffi- 
cult to observe. Hence a model which neglects screening, 
like the gauge glass, should be able to account for most 
of the observable data in the critical region. 

Unlike a real superconductor in a magnetic field, the 
gauge glass is isotropic on average. There are local 
quenched fluxes but no net field in any direction. In six 
or more spatial dimensions the lack of anisotropy does 
not seem to matters; however it is still an open question 
whether, in three dimensions, the vortex glass transition 
in the gauge glass and in a system with a non-zero net 
field are in the same universality class. 

The source of the quenched randomness in the gauge 
glass model is the vector potentials linking the sites. This 



is not very realistic, and a more accurate model would 
have vector potentials corresponding to a uniform field, 
and put disorder into the strength of the couplings. How- 
ever, the precise details of the disorder should be irrele- 
vant for critical phenomena 

In this paper we shall show very clearly that the gauge 
glass has a finite temperature transition to a vortex glass 
ordered state and that its critical exponents agree, within 
fairly large error bars, with some experimental measure- 
ments on type-II superconductors. More work remains 
to be done to check whether an anisotropic model with 
a net field would change the universality class. 



III. THE OBSERVABLES 



A standard technique to detepaaine the critical temper- 
ature is to use the Binder raticO 



g(L) = 2- 



[<M 4 >] a 



M 2 >] 



(3) 



where [• ■ -] av denotes an average over configurations of 
the disorder, (• • •) denotes a thermal average and q is the 
complex overlap order parameter 



1 



N 

V. 



1 = ^2^ ex PW 

3=1 



(4) 



in which a and (3 denote two independent replicas with 
the same disorder. One plots g(L) vs. T for different 
L and identifies the temperature at which the curves 
cross as T c . However, since g{L) cannot exceed unity, 
the splaying out of the data for g(L) below T c (which 
indicates spin glass order) is a small effect which can be 
difficult to see. This is why it as been so hard to establish 
conclusively that there is spin glass c*dej^,at finite-T in 
the three-dimensional Ising spin-glass Jl3il3 ._. 

In order to avoid this problem we follow Reger et al.E3 
in calculating a current; this is defined as the rate of 
change of the free energy with respect to a twist angle 
O at the boundaries. We begin by replacing periodic 
boundary conditions with twisted boundary conditions 
along one of the axes x, i.e. 



e. 



(5) 



Note that 9 = corresponds to periodic boundary condi- 
tions, and O = 7r corresponds to anti-periodic boundary 
conditions. We can convert the model back to periodic 
boundary conditions by redefining <fi through 



J i+nx 



J i-{-nx 



(6) 



The model then is precisely Eq. ([!]) with periodic bound- 
ary conditions but with the Aij for bonds in the x- 
direction changed according to 



2 



e 

-<*-i,i-\-x j • 



(7) 



Using F = — ln(Z) / f3 we define a current as the response 
of the free energy to an infinitesimal 0, 



OF 

I(L) = lim 



.4 



Because the A^- are uniformly distributed over the en- 
tire period of the sine function, the the value of the cur- 
rent averaged over samples is zero, i.e. 



[I(L)] av = 0. 



(9) 



Consequently we calculate the root-mean-square current, 
given by 



y/[{I a (L))(Ip(L))] a 



(10) 



where a and /3 denote two independent replicas with the 
same disorder. We use two replicas to avoid any bias in 
calculating the average of the square of the current. 

The primary advantage of using 7 rms rather than g(L) 
is that I rms increases with L for T < T c , so the splaying 
out the data below T c should be much easier to see than 
for the Binder ratio. In the ordered phase the current 
should scale with the stiffness exponentt3E-at3 9, i.e. 



in 



(T < T c ), 



(11) 



where 9 > if T c > 0. Above T c , where spin-spin corre- 
lations are small, we expect J rms to approach zero with 
increasing L because larger systems are less sensitive to 
a twist at the boundaries. The J rms curves cross at T c 
and have the finite size scaling form 



where 



t=(T— T c )/T c 
is the reduced temperature. 



(12) 



(13) 



EXCHANGE MONTE CARLO AND 
EQUILIBRATION 



IV. 



Exchange-, Monte Carlo,ElE!J also called parallel 
tempering,c2Tc3 is a technique for simultaneously simu- 
lating multiple copies of a particular configuration of dis- 
order with each copy at a different temperature. After a 
certain number of sweeps through the lattice, one tries to 
exchange the spin configurations of copies at neighboring 
temperatures with the probability 



P(cr m <-» <7 m+ i;/3 m ,/3 m+ i) = exp(-A) 



(14) 



{fi m - p m+1 )(E(a m+1 ) - E(a m )). 



(15) 



<Jm is the spin configuration of the mth copy which has 
inverse temperature (3 m and total energy E(<r m ). A given 
spin configuration is thus heated and cooled many times 
during the simulation. Since equilibration is fast at high 
temperature, each time the system is cooled the mini- 
mum (valley) that it visits is uncorrelated with the min- 
imum that it visited at the previous cooling. Hence the 
system can can visit different local minima at low temper- 
ature more efficiently than if the temperature were kept 
fixed. In the latter case, very large barriers would have 
to be overcome, which takes a time exponentially large in 
the ratio of the barrier height to the temperature. 

In deciding what temperatures to simulate, one would 
like the energy distributions of neighboring temperatures 
to have enough overlap that the probability to exchange 
configurations is sufficiently high. This requires 



T, 



m+l 



T 



(16) 



where the c m are constants of order unity. In our sim- 
ulations we took all the c m to be exactly unity, which 
gave a satisfactory acceptance ratio for temperature ex- 
changes in the interval from 0.5 to 0.6 for all sizes and 
temperatures studied. 

We took the highest temperature, T max , to be approx- 
imately 2T C at which the spin dynamics are quite rapid. 
Temperature exchanges were carried out after every 10 
Monte Carlo sweeps, after which "time" the normalized 
energy-energy auto-correlation function at T max is quite 
low, about 0.25. 

The equilibration time, t oq , for I ims is determined by 
following the temporal evolution of 



I 2 (t ) 



^ to 

-^I a {t + t )Ip{t + t ) 



to 



4=1 



(17) 



where to is the number of equilibration sweeps as well the 
number of measurement sweeps and the subscripts a and 
P denote independent replicas. When to <C t Qq the spin 
configurations of the two replicas are completely uncor- 
related but as to increases they become more correlated, 
since they both feel the same random interactions. Thus 
we expect I 2 (to) to monotonically increase from zero to 
the equilibrium value as to — » ^c q ; see Fig. (||). Each of 
the equilibration times in Table | is chosen to be the least 
number of sweeps necessary to equilibrate at the lowest 
temperature, T m ; n . For each size except L = 12 we ran 
some samples for to 3> t eq in order to confidently deter- 
mine t eq ; the remaining samples were run with to = t oq . 
We should note that for L = 12 it was impractical to 
run jobs for to much greater than the value at which the 
J 2 (to) seemed to stop increasing, so we took this value to 
be t eq . 



where 
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T ■ 

' min 


T 

- 1 max 




t 




4 


0.3 


0.97 


n 


1280 


8000 


6 


0.3 


0.92 


18 


10240 


10937 


8 


0.3 


0.92 


27 


20480 


5388 


12 


0.3 


0.90 


47 


163840 


781 



TABLE I. Parameters of the exchange Monte Carlo sim- 
ulations for each value of L, where T m i n and T max are the 
minimum and maximum of m temperatures, and t eq is the 
number of sweeps for equilibration, which, in our simulations, 
is also equal to the number of sweeps for measurements. The 
number of samples studied for each size is also shown. 



which agreeg-with the previous value of 0.45 ± 0.05 from 
Reger et al.E^I The key difference is that in Reger et al. 
the values for J rms did not splay out significantly below 
T c and so they did not find compelling evidence for spin 
glass order. 



L = 4 


L = 6 


L = 8 


L = 12 


1.026 ± 0.009 


1.090 ± 0.008 


1.165 ±0.012 


1.238 ± 0.037 



TABLE II. Jrms from Eq. (JlC|) for T/J = 0.3 and differ- 
ent L. The values are clearly distinct beyond the size of the 
error-bars. 
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FIG. 1. A semi-log plot of I (to) against to, from Eq. (|17 
for L = 4atT/J = 0.3 using 2000 samples to test for equili- 
bration. One sees that the results do not change for to greater 
than about 1000. Hence, in the production runs, the equili- 
bration time was taken to be to = 1280. 
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FIG. 2. A plot of J rms from Eq. (|10j) for L — 4, 6, 8, and 12. 
For the point of intersection we estimate that the critical tem- 
perature is 0.47 ± 0.03. Some of the L — 12 errorbars have 
been removed for clarity. 



V. RESULTS FOR STATICS 

We present our results for 7 rms from Eq. ([l0]) in 
Fig. (g). The data cross at T c /J — 0.47 and splay out 
clearly on both sides of the transition. This is the first 
time that the splaying out of the data for I rms in the or- 
dered state can be distinguished well beyond the size of 
the error bars, thus presenting incontrovertible evidence 
for a spin-glass ordered phase in three dimensions. The 
data for the lowest temperature are given numerically in 
Table (||). We obtain an uncertainty of 0.03 for T c /J by 
setting the error equal to the region over which the data 
for all sizes overlap within their error bars. Our final 
estimate for the critical temperature is 



T c /J = 0.47 ±0.03, 



(18) 



We obtain an estimate for the correlation length expo- 
nent v by finite-size scaling the data from Fig. (Q); see 
also Eq. (|l2|). Fig. (||) shows that the data scales well 
with v = 1.39. This value is obtained by by scaling the 
data using different choices for z/, fitting a polynomial, 
calculating the chi-squared statistic for each choice, and 
minimizing the chi-squared. The error in v is determined 
by varying the value of v until the L — 8 data no longer 
scales with the smaller sizes within the sum of their er- 
rorbars. L = 8 was chosen rather than L = 12 because 
the errorbars of the former are much smaller. This leads 
to our estimate 



v = 1.39 ±0.20. 



(19) 



Variations in T c from Eq. (18) did not significantly in- 
crease the error bar beyond what that shown in Eq. (|i"s|). 
Our result agrees with, and is a bit more precise than, 



4 



those previously presented bv-Wengel and YoungJ^ v 
1.3 ± 0.3, and by Reger et al,E3 1.3 ± 0.4. 

The correlation length exponent has also been deduce- 
from experimental measurements. It is now knownjj' 
however, that the response of the resistivity and crit- 
ical temperature to tilting of the applied field distin- 
guishes a vortex-glass transition, in which point disorder 
dominates, from a Bose-glass transition, in which cor- 
related disorder, such as twin boundaries or columnar 
pins, dominates. To our knowledge, there are only two 
experimentaj'EH which observe a vortex glass transition 
and demonstrate the proper response to magnetic field 
tilting. p 

In the first, Petrean et alJ9 find that for untwinned 
proton-irradiated YBa2Cu307_5 near criticality the re- 
sistivity decreases from its maximum as the field is tilted 
away from the c axis,El, which signals a vortex glass tran- 
sition, whereas in twinned YBa2Cu307_5, they find a 
Bose-glass transition with thejjesistivity increasing from 
its minimum.Q Petrean et al.Jj do not obtain v but do 
obtain the resistivity exponent s by fitting resistivity vs. 
temperature curves. We will compare our results to these 

In the second work, Klein et al.,EZI study the vortex 
glass transition in (K, Ba)BiC>3. The dependence of their 
data on field tilting is that expected for the vortex glass, 
and they obtain v = 1.0 ± 0.2 which agrees with ours 
within the erparbars. 

KawamuraEj has recently modified the gauge glass to 
include a net magnetic field. He obtains v — 2.2 ± 0.4 
which is greater than, our estimate and the experimental 
result of Klein et a£j. 




FIG. 3. A plot of / rms from Eq. (|T(J) versus L 1/v t, where 
t = (T - T c )/T c and v = 1.39. Some of the L — 12 errorbars 
have been removed for clarity. 



In addition to v, we calculate the independent expo- 
nent r\ which describes the decay of the correlation func- 
tion at criticality. To our knowledge, this is the first 
calculation of r\ for the gauge glass in three dimensions. 
Finite size scaling predicts that, at the critical point, the 
spin glass susceptibility 



XSG = N[(\ q \ 2 )} : 



scales as 



XSG 



L 



2-77 



(20) 



(21) 



assuming the hyperscaling relation j/v = 2 — rj. We 
use standard Monte Carlo to calculate the susceptibility 
at T/J = 0.44, 0.47, and 0.50; the details of these simu- 
lations are discussed in § VI . We present our results in 



Fig. (|4|). From linear least-squares fits to the data on a 
log-log scale we obtain 



rj = -0.47 ±0.07. 



(22) 



The error in our estimate-comes mainly from the uncer- 
tainty in T c . Kawamura'sEj anisotropic gauge glass yields 
a similar estimate, rj = —0.5 ± 0.2. 



o 
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10 



L 



FIG. 4. A log-log plot of xsg vs. L for L = 4, 6, 8, and 10 
at T/J = 0.47 according to Eq. @ with 77 = -0.47. 

There has been some controversy regarding the value 
of the zero temperature stiffness exponent 9. This is typ- 
ically computed from the root-mean-square "defect en- 
ergy" on changing the boundary conditions from periodic 
to anti-periodic, but is also given by the response to an in- 
finitesimal twist in the boundary conditions, as shown in 
Eq. dll|). A positive value for 9 indicates a finite tempera- 
ture transition to a spi&|g4|asSr*tate, whereas 9 < implies 
T r = 0. Some studiesNI^By obtain < 9 < 0.077; we 



5 



refer to these as the low group. Others however, 
calculate 0.26 < 9 < 0.31; we name these the high group. 

Our results for the rms current are not at sufficiently 
low T and large L to get a firm estimate for 9, but we 
obtain a plausible bound as follows. For each tempera- 
ture below T c we do a linear least squares fit of log 7 rms 
against logL to get an effective temperature dependent 
value, 9 e ff(T). The results are shown in Fig. (||). We 
see that 9 e s increases monotonically as T decreases, so 
we expect that the asymptotic value is greater than the 
value of 9 e s at the lowest temperature, i.e. we expect 

(9 > 0.18, (23) 

which is consistent with the results of the "high group" . 
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< 0.1 
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0.3 0.35 0.4 0.45 
T/J 

FIG. 5. A plot of e cff (T) vs. T. At each T, 9 e g(T) 
is obtained from a linear least squares fit of log/ rms against 
logL. 



VI. DYNAMICS 

When the gauge glass is sufficiently close to criticality 
the correlation length is of order L and the system expe- 
riences critical slowing down. The equilibration time t cq 
then scales as 



eq 



L\ (24) 

where z is the dynamic exponent. At temperatures just 
above the vortex glass transition, the resistivity p is 
predictedQ to scale as 

p~(T-T c ) s , (25) 

where s is related to other exponents by 



v(z-l) 



(26) 



in three dimensions. We are thus motivated to calculate 
z to compare our results to experimental measurements 
of s. It must be emphasized, however, that there are 
more dynamical universality classes than static univer- 
sality classes. In our simulations which determine z we 
use dissipative (Monte Carlo) dynamics with standard 
Metropolis-type updating probabilities, and without the 
temperature swapping (exchange Monte Carlo) that we 
used in our simulations of static quantities. Tests lor 
equilibration were carried out as in Bhatt and YoungtLj. 



T/J 


L = 4 


L = 6 


L = 8 


L = 10 


0.44 


500 


400 


358 


372 


0.47 


500 


400 


328 


246 


0.50 


500 


400 


350 


321 



TABLE III. Number of samples in the standard Monte 
Carlo simulations used for calculating the dynamic exponent. 
Three temperatures were used to account for the uncertainty 
in T c . 

The dynamic exponent can be obtained by a finite size 
scaling of time-dependent measurements of the spin glass 
susceptibility. We use the "two-replica" susceptibility 
Xsg(£o) defined by 



Xsg(*o) 



N 

y^ e «[^(to+t)-^(*o+*)] 

i=i 



(27) 



where a and (3 denote independent replicas. 

We obtain our estimate of z by a finite-size scaling 
analysis of the data for xsg(^o) for sizes L = 6, 8, and 10; 
the scaling is better without the L = 4 data, so it has 
been omitted. One expects from Eq. ( |2^ ) that at criti- 
cality 



Xscfa) 



f(L~ z t ) 



(28) 



cq ) 



where / is a scaling function. Our results, shown in 
Fig. (|), yield 



4.2 ±0.6, 



(29) 



as calculated from the T/J = 0.47 data. The estimate of 
the uncertainty is obtained by varying z until the L = 10 
data no longer scales with the smaller sizes within the 
errorbars. Similar scaling plots from the T/J = 0.44 
and T/J = 0.50 yield z values within Eq. (|29|), thus the 
uncertainty in T c is not the dominant contribution to the 
error in z. Our value agrees path the previous estimate 
z = 4.7 ± 0.X|by Reger et alE3 and z — 3.3 ± 0.5 from 
Kawamura'sH anisotrop ic g auge glass. 

Using Eqs. ( p9|) and (|i~9|), the resistivity exponent is 
obtained from Eq. 



(26) 
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FIG. 6. A finite-size scaling plot of xsG(to) using Eq. (|2S 
for L = 6,8, and 10 at T/J = 0.47 with z = 4.2. 



4.5 ±1.1 



(30) 



This value agrees with the experiments of Klein et al.l 
s = 3.9 ± 0.3, on (K, Ba)Bi03, as well as Petrean et al. 
s = 5.3 ± 0.7, who study untwinned proton-irradiated 
YBa 2 Cu 3 7 _ 5 . 



VII. CONCLUSIONS 

In this work we have presented exchange Monte Carlo 
results for the gauge glass model in three dimensions. We 
have shown for the first time incontrovertible evidence of 
a vortex-glass ordered phase below T c . To our knowl- 
edge, this is the first time that such clear ordering has 
been shown for any spin-glass model in three dimensions. 
The correlation length exponent v has been calculated to 
higher precision than before. We have also obtained the 
first estimate of the correlation function exponent rj from 
Monte Carlo simulations. Finally, our values of v and z 
are combined to estimate the resistivity exponent s. Our 
results are summarized in Table ( [V ) . 

It is interesting that our values of rj and v agree 
with those of two other glassy systems, the three- 
dimensional ± J Ising spin glassliil and chical ordering of 
the three-dimensional ±JXY spin glassH Kawashima 
and YoungM obtain v = 1.7 ± 0.3 and n rn-0.35 ± 0.05 
for the Ising spin glass, while KawamuraES obtains v = 
1.5 ± 0.3 and 77 = —0.4 ± 0.2 for the chiral glass order- 
ing. However, it is expected that the Ising and gauge 
glass models are in different universality classes because 
the order parameter has a different number of compo- 
nents in each case: one for the Ising spin glass, and two 



(the real and imaginary parts of Eq. (H)) for the gauge 
glass. Indeed, a first order e — 6 — d expansion shows 
that the gauge glass exponents-are not the same as any 
n-component vector spin glass.Ej The order parameters 
for the Ising spin glass and the chiral transition in the XY 
spin glass do have the same dimensionality, but the tran- 
sitions are still likely to be in different universality classes 
because of the long range interaction between chiralities 
in the XY spin glass. The error bars in the exponents are 
not extremely small so the apparent agreement between 
the results may be simply a numerical coincidence. 



T c /J 


V 


V 


z 


0.47 ±0.03 


1.39 ±0.20 


-0.47 ±0.07 


4.2 ±0.6 



TABLE IV. Critical temperature and exponents of the 
gauge glass in three dimensions. 

To confidently compare these results with experiment 
one has to show that the critical point in the gauge glass 
model is in the same universality class as an anisotropic 
model with a net field. Clearly, if such anisotropy causes 
the scaling behavior to be anisotropic, in the sense that 
the exponents for the divergence of the correlation length 
are different for separations along and perpendicular to 
the field, then it is relevant. If, however, the anisotropy 
induced by a net field does not lead to anisotropic scaling, 
but just causes the amplitudes of the correlation lengths 
in the different directions to be different, then it is not 
clear why it should be relevant, just as making an Ising 
ferromagnet anisotropic by having the bonds in one di- 
rection different from those in the other directions does 
not change the. universality class. 

KawamuraEJ has studied a vortex glass model with a 
net field and finds different critical behavior from that of 
the gauge glass, even though the scaling is not anisotropic 
in the sense defined above. However, he imposes free 
boundary conditions which may lead to very large cor- 
rections to finite-size scaling (since a large fraction of the 
sites are on the boundary). Hence it would be very useful 
to study the vortex glass transition in a model with a net 
field and periodic boundary conditions. 
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